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ABSTRACT:    Area-elevation  information  can  be 
collected  -with  grid  intersection  sampling  techniques 
more  rapidly  than  by  planimetering.    A  nomograph 
speeds  the  selection  of  grid  spacing  for  l)  various 
sampling  intensities,  2)  drainage  areas  from  1  - 
10,000  square  miles,  and  3)  five  standard  topographic 
map  scales.    Routine  calculations  and  linear  inter- 
polations within  elevation  classes  are  carried  out  by 
an  IBM  1620  computer  using  a  FORTRAN  program. 


Area-elevation  curves  have  proved  useful  in  evaluating  a  variety 
of  meteorological  and  hydro logic  conditions  (Anderson  1958 ;  Corps  of 
Engineers  1956;  Langbein  19^7;  Strahler  1952;  Soil  Conservation  Service 
1957)-     These  uses  range  from  temperature  extrapolations  to  the  determina- 
tion of  reservoir  storage  capacities.     One  author  has  even  used  such  curves 
to  study  fluctuations  in  glacial  regimen  (Baussart  1961). 

Planimetering  and  grid  sampling  are  the  two  accepted  data-collec- 
tion methods  for  establishing  area -elevation  relationships.  Planimeter- 
ing the  area  between  contour  lines  is  the  standard  and  most  commonly  used 
technique  for  Federal  and  private  water  resource  construction  agencies. 
It  provides  a  complete  survey  of  the  watershed  area  but  is  slow  and  sub- 
ject to  human  and  instrument  errors  in  the  measurements. 

In  contrast,  grid  intersection  sampling  is  faster  and  easier  to 
apply  than  planimetering.     Linsley,  Kohler,  and  Paulhus  (19^9)  state 
that  it  yields  "equally  satisfactory  results"  and  that  "In  general,  at 


least  100  points  should  be  used. "    Still  the  method  has  not  gained  much 
favor  with  hydrologists  in  the  30  years  since  it  was  proposed  by  Horton 
(1932).    Perhaps  skepticism  about  the  accuracy  of  sampling  explains  this 
lack  of  acceptance. 

Our  preliminary  tests  have  shown  that  area -elevation  curves  con- 
structed by  grid  sampling  will  coincide  with  planimeter -determined  curves , 
if  enough  sampling  points  are  taken.     Usually  less  than  500  points  are 
needed  with  8  elevation  classes.    As  few  as  100  points  will  yield  highly 
useful  information  for  most  purposes.     The  two  techniques  are  compared  for 
the  South  Fork  Eel  River  above  Miranda,  California  (table  l).     Grid  sampling 
results  are  based  upon  108  intersections.     The  maximum  difference  in  percent- 
area  within  a  class  was  3 • 57  (19  square  miles).     The  difference  between  esti- 
mated mean  basin  elevations  was  85  feet. 

For  many  uses  these  differences  would  not  be  critical.     For  example, 
temperature  extrapolations  to  the  mean  basin  elevation  would  differ  only 
0.25  to  0.35°F.,  depending  upon  what  lapse  rate  is  assumed.     If  the  limits 
of  permissible  error  are  more  stringent,  greater  accuracy  may  be  obtained 
by  increasing  the  sampling  intensity. 

Table  1. --Area -elevation  table,  South  Fork  Eel  River  above  Miranda, 
California — area,  537  square  miles 


Elevation  class, 

Perc ent 

area  within 

elevation  class 

feet 

\  Grid 

Planimeter 

Difference 

250-1,000 

20.37 

18.90 

1.^5 

1,000-1,500 

23.15 

20.09 

3.06 

1,500-2,000 

23°  15 

25.07 

1.92 

2,000-3,000 

28.70 

27.79 

•91 

3, 000 -U, 491 

k.63 

8.20 

3-57 

Mean  elevation: 

...1,798  ft. 

Grid   0 

...1,713  ft. 

In  grid  sampling,  a  uniform  grid  is  laid  over  a  topographic  map(s) 
of  the  watershed.     Considerable  time  can  be  expended  determining  the  proper 
spacing  of  the  grid  lines  for  a  given  sampling  intensity  as  this  is  a  func- 
tion of  both  drainage  area  and  map  scale.    A  nomograph  (fig.  l)  eliminatei- 
this  bottleneck.    To  use  the  nomograph,  enter  with  the  drainage  area  of  the 
watershed  (a)  and  move  vertically  until  the  intersection  with  the  sampling 
intensity  line  (n).     Then  move  horizontally  to  the  appropriate  scale  line 
(S)  and  vertically  to  the  grid-spacing  scale  (i).    For  example,  we  have 
determined  the  grid  spacing  for  the  South  Fork  Eel  River  with  100  sampling 
points.    The  nomograph  shows  that  the  spacing  should  be  2.37  inches  (scale  l)« 
This  would  be  a  cumbersome  interval  so  every  other  line  on  a  1  inch  grid  was 
used. 
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Between  7  and  1^  elevation  classes  are  established  on  a  watershed. 
The  exact  limits  of  each  class  are  determined  by  the  contour  interval  and 
the  technician's  judgment  of  the  elevation  distribution  within  the  watershed. 
Ideally,  about  the  same  number  of  points  should  fall  within  each  class. 
Colored  pencils  can  be  used  to  trace  class  boundaries  to  aid  in  distinguish- 
ing the  different  classes.    Each  grid  intersection  is  then  tallied  within 
the  appropriate  elevation  class. 

All  subsequent  calculations  are  made  on  an  IBM  1620  computer  with 
a  20K  memory.    A  FORTRAN  program  (fig.  5)  for  the  computer  is  general  for 
7  to  Ik  elevation  classes  and  for  any  type  of  relative  area  measurement 
(either  sampled  or  planimetered) . 

Input  for  the  program  (figs.  2  and  3)  consists  of:     (l)  a  watershed 
identification,  (2)  drainage  area,  (3)  maximum  basin  elevation,  (k)  mini- 
mum elevation  of  each  class,  and  (5)  an  integer  value  of  relative  area 
within  each  class.     The  watershed  may  be  identified  by  any  series  of 
thirteen  alphanumeric  characters.     For  example,  if  the  U.S.  Geological 
Survey  gage  number  is  used,  three  digits    will  remain  for  dummy  numbers 
or  regional  codes.     The  drainage  area  may  range  from  .000001  to  99,999-99 
square  miles. 

Output  (fig.  k)  includes:     (l)  class  range,  (2)  class  mean,  (3) 
percent  area  in  each  class,  (k)  mean  watershed  elevation,  and  (5)  an  area- 
elevation  table.     The  area-elevation  table  lists  elevations  by  0.1  increments 
of  the  original  elevation  classes,  with  the  cumulative  area  in  square  miles, 
and  cumulative  percent  area  below  that  elevation.     Linear  interpolation  within 
each  class  is  used  because  it  is  the  simplest  model  and  there  is  no  theoret- 
ical justification  for  any  curvilinear  function  within  a  class.    An  output 
example  is  appended. 

In  practice  we  have  usually  required  that  the  lowest  and  highest 
classes  have  at  least  one  sample  point.     If  no  sample  points  fall  within 
the  highest  or  lowest  class,  we  combine  them  with  the  next  adjacent  class. 
This  procedure  avoids  nonsense  statements  in  the  table,  such  as  "100.00 
percent  of  the  area  is  below  9; 000  feet"  when  the  highest  basin  elevation 
is  actually  9,500  feet. 

If  a  1620  computer  is  available  with  a  larger  memory  capacity  than 
20K,     the  program  could  be  modified  so  that  a  larger  number  of  elevation 
classes  could  be  used.     For  most  problems  additional  classes  are  not 
necessary.    A  second  program  modification  would  be  warranted  if  a  1620 
Auto -Plotter  were  available.     This  attachment  could  be  programed  to 
plot  the  tabular  material  of  the  area-evaluation  calculations.     In  this 
case  also,  a  larger  memory  unit  would  be  necessary. 
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Figure  2. --Data  input  format. 


Columns  Fields  Title 

Data  Header  Card 

1-k  (l^-)  Number  of  watersheds  to  be  processed. 

Data  Card  One 

h- 16  (l3H)  Alphanumeric  watershed  name. 

17-29  (F13-6)  Watershed  area  in  square  miles. 

32-37  (F6.0)  Maximum  basin  elevation. 

k'i-kk  (12)  Number  of  elevation  classes. 

Data  Card  Two  (One  per  watershed) 

19-23  Cf^.O)  Minimum  elevation  of  1st  class. 

27-31  (f5»0)  Minimum  elevation  of  2nd  class. 

35-39  (F5«0)  Minimum  elevation  of  3rd  class. 

^3-^7  (F5-0)  Minimum  elevation  of  kth  class. 

51-55  (F5«0)  Minimum  elevation  of  5th  class. 

59-63  (F5«0)  Minimum  elevation  of  6th  class. 

67-71  (F5.0)  Minimum  elevation  of  7th  class. 

Data  Card  Three  (One  per  watershed) 

19-23  (F5-0)  Relative  area  in  1st  class. 

27-31  (F5.0)  Relative  area  in  2nd  class. 

35-39  (F5.0)  Relative  area  in  3rd  class. 

^3-^7  (F5.0J  Relative  area  in  kth  class. 

51-55  (F5'0)  Relative  area  in  5"th  class. 

59-63  (F5-0)  Relative  area  in  6th  class. 

67-71  (F5-0)  Relative  area  in  7th  class. 


Data  Card  Four  (One  or  none  per  watershed) 

Same  format  as  Data  Card  Two  but  minimum  elevations  are  for 
elevation  classes  8  to  Ik.     If  less  than  1^  classes  leave  remaining 
fields  blank. 

Data  Card  Five  (One  or  none  per  watershed) 

Same  format  as  Data  Card  Three  but  relative  areas  are  for  classes 
8  to  1^.     If  less  than  Ik  classes  are  used  leave  che  remaining  fields 
blank. 
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Figure  3. 

EXAMPLE  OF   AREA-ELEVATION  DATA  INPUT 


DATA  CARD 
DATA  CARD 
DATA  CARD 
DATA  CARD 
DATA  CARD 
DATA  CARD 
DATA  CARD 
DATA  CARD 


HEADER  CARD 

F91220 

F83100 


F91600 


12.5 

335 
24 
105. 
50 
T7— 

2800 

_  r 

82. 

380 
2 


2750 
700  1000 
11  18 
3209 
400  800 
TT  28 


800 
7 


4480 
1300 
11 


7 

1300 
7 

8 

1200 

~20 


1800 
13 


1600 

8 

1600 


2300 
20 


1900 

6 

2000 
3 


2800 
12 


2200 
1 

2400 

~~ 2 


3  30,0 
3 


  Figure  k.   

EXAMPLE  OF   DATA  OUTPUT  FORMAT 


WATERSHED 
NAME 

F91220 


AREA  IN 
SQU  MI 

12.500000 


MEAN  BASIN 
ELEV   IN  FEET 

1085. 


MIN     CLASS  ELEV  335.  700.  1000.  1300. 

MAX     CLASS  ELEV  700.  1000.  1300.  1600. 

MEAN  CLASS  ELEV  518.  850.  1150.  1450. 

PERCENT   OF   AREA  32.00  14.67  24.00  9.33 


1600.  1900.  2200. 

1900.  2200.  2750. 

1750.  2050.  2475. 

10.67  8.00  1.33 


CUMULATIVE  AREA  AND  CUMULATIVE  PERCENT 
AREA  BELOW  LISTED  ELEVATIONS 

ELEVATION  CUMULATIVE  CUMULATIVE 

IN   FEET  SQU  MILES  PERCENT 

371.  .40000  3.20 

408.  .80000  6.40 

444.  1.20000  9.60 


2640.  12.46665  99.73 

2695.  12.48332  99.86 

27507"  T2. 49998"  99. 99~ 
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Figure  5. 

FORTRAN   SOURCE   PROGRAM  FOR  1620 


1 

r ORM A  I 

1  DA  » 

34HWATERSHED 

ADC  A       T  M 

AKtA     1  IN 

R  A  c.  I  M  \ 
D  Mo  1  IN  J 

2 

r™  /~\  nil  A  T 

FORMAT 

(  9X  ♦ 

33HNAME 

SUU    M 1 

PI    P  W  HI 

t  L  L  V  IN 

P  P  P  T  \ 

r  1 1  1  1 

3 

r ORM A  T 

(IX) 

5 

FORMAT 

(  1  6H 

MI  N 

CLASS 

ELEV7F8.0 ) 

6 

FORMAT 

(  16H 

MAX 

CLASS 

ELEV7F8.0 ) 

7 

FORMAT 

(  1  6H 

MEAN 

CLASS 

ELEV7F8 .0 ) 

9 

FORMAT 

(  1  6H 

PERCENT  OF 

AREA7F8  •  2  ) 

10 

FORMAT 

(  43H 

CUMULATIVE   AREA  AND 

CUMULATIVE  PERCENT  ) 

1 1 

1 —  a  n  ii  a  t 

FORMAT 

/   O  /  L_J 

(  34H 

AREA   BELOW  LISTED 

ELEVAT IONS ) 

1  2 

FORMAT 

(  43  H 

ELEVATION  CUMULATIVE 

rriiurrn  a  t  t  \/c  i 

LUMUL A  1  1 Vt ) 

13 

FORMAT 

(  40H 

T  M  ccct 

1  in    r  t  L  1 

SQU  MILES 

PERCENT ) 

14 

FORMAT 

(  6X  ♦ 

r  /  •  0  » 

0  A  »  r  1  £  • 

_>*M-A»ro»l) 

1  5 

r ORM AT 

(  Ion 

M  i  IN 

ri  Ace 

ELEV7F8 • 0 ) 

16 

FORMAT 

/    *1    /  LJ 

(  1  6H 

M  I  N 

CLASS 

tLLVoroiU 1 

1 7 

FORMAT 

11/11 

(  16H 

M  T  M 
M  1  l\ 

C  1   A  C  c 

pi  pwcpo .n l 

1  fl 

J.  O 

FORMAT 

HAH 

M  T  M 
H  1  IN 

C  \  Ace 

pi  ri//.  pa    a  1 

1 9 

r-  a-\  r~\  ft  i  A  T 

r ORM AT 

(  16H 

M  1  IN 

C I  A  c  c 

ELEV3F8 • 0 ) 

20 

r  Ah  y  A  T 

r ORMAT 

(  1  6H 

M  I  N 

CLASS 

PI  pi/t  PD    n  1 

2  1 

r ORMAT 

(  1  6  H 

MIN 

CLASS 

ELEVF8 . 0 ) 

22 

FORMAT 

(  16H 

MAX 

CLASS 

PI  F\/7Pfl     n  | 
tLtv  1  ro • U 1 

£  3 

r UKM A  1 

I  1  O  H 

MAX 

CLASS 

ELEV6F8 . 0 ) 

"3  /• 

r UKM A  1 

lion 

MAX 

CLASS" 

pi    CWC  CD  ~  f\  \ 

1  UKMM  1 

lion 

MAX 

CLASS 

pi  P\//iPfl.nl 

26 

FORMAT 

{  16H 

MAX 

CLASS 

cri  cx/qPd  n\ 

9  "7 

rUKrlH  1 

lion 

MAX 

CLASS 

ELEV2F8.0) 

~)  D 
£  O 

lion 

MAX 

CLASS 

Pi  P \ / P  c    n  \ 
tLtVi  0  »(j  ) 

^  9 

CT  f*l  D  KA  A  T 

r UKMA 1 

I  I  6H 

MEAN 

CLASS 

ELEV7F8 • 0 ) 

'   1 —  1 —  v     1    1     WW  \J  f 

3n 

FORMAT 

(1  6H 

MEAN 

CLASS 

PI   P  \  /  c  P  0     r\  i 

ELt V6F  0 • 0 ) 

3  1 

r  Ar-)M  A  T 

F ORMAT 

(  1  6H 

MEAN 

CLASS 

ELEV5F8 .0 ) 

3? 

FORMAT 

{ i  fsH 

V   J>  u  1  ' 

MEAN 

CLASS 

ELEV4F8.0 ) 

r ORMAT 

/   1   /*  LJ 

(  1 6H 

MEAN 

CLASS 

ELEV3F8 .0 ) 

34 

FORMAT 

(  1  6H 

MEAN 

CLASS 

ELEV2F8.0) 

r UKMA I 

lion 

MEAN 

CLASS 

ELEVF8.0 ) 

36 

F ORMAT 

(  1  6H 

PERCENT  OF 

AREA7F8.2 ) 

Q  "7 

3  7 

r ORMAT 

#   1   r  \  1 

(  16H 

PERCENT  OF 

AREA6F8.2 ) 

3ft 

FORMAT 

PERCENT  OF 

AREA5F8 • 2 ) 

cr  r\  r>  u  a  x 

FORMAT 

(  16H 

PERCENT  OF 

AREA4F8.2 ) 

FORMAT 

(  1  6H 

PERCENT  OF 

AREA3F8.2 ) 

4  1 

FORMAT 

(  16H 

PERCENT  OF 

AREA2F8.2 ) 

42 

P  r\  hi  KJ  A  T 

FORMAT 

(  16H 

PERCENT  OF 

ARE AF8 • 2  ) 

43 

FORMAT 

(  14  ) 

44 

FORMAT 

(  3X  ♦ 

13H 

NAME,F13.6» 

2X»F6.0» 

5X» 12 ) 

45 

FORMAT 

(  15X 

♦  7F8.0  ) 

46 

FORMAT 

(  6X  » 

F6.  0  ♦ 

5X.F12. 

5  »5XtF7.2  ) 

47 

FORMAT 

(15HMISTAKE    IN  DATA) 

48   FORMAT    (43HNEXT   DATA  CARD  MUST   BE  FORMAT   44  *  PUSH  START ) 


DIMENSION  E(15)»R(15)»P(14),S(14),D(1A),F(14) 

NOS  =  0 
J  =  l 

READ  43.NUM 
KACC=0 
89  PAUSE 
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Figure  5. —Continued 


100  READ  44»AREA»EMAX»JAC 

KACC=KACC+1 

DO   90  J=ltl5 

90  R( J)=0. 
J=l 

IF    ( JAC-14)50t50»114 
50    IF    (JAC-7)  114,102,101 

102  READ  45  ,E(  1  )  ,E( 2 }  ,E(3 J ,E(4) ,E(5 ) ,E(6 ) ,E(7) 
READ  45,R(1),R(2),R(3),R(4),R(5),R(6),R(7) 

99  E{ JAC+1 )=EMAX 
RCON=R ( JAC+1 ) 
IF   (RCON) 114,91 .114 

91  R AEA=0 • 
SUME=0. 

DO   103  J=1,JAC 

103  RAEA=RAEA+R( J) 
DO   98  J=1,JAC 

P( J)=(R( JJ/RAEA)  

D( J)=E( J+1)-E(J) 
F(J)=D(J)/2.+E(J) 
SUME=SUME+(F( J)*P( J) ) 
P( J)=P( J J*i00.+.005 
IF    (D( J) )        114» 114,98 

98  FU)=F(J>  +  .5 —  

PUNCH  3 
PUNCH  3 
PUNCH  3 
PUNCH  1 
 PUNCH  2 

PUNCFT3  "~  

PUNCH  44»AREA,SUME ,NOS 
PUNCH  3 
PUNCH  3 

PUNCH  5»E(1),E(2)»E(3),E(4J.E(5),E(6),E{7 ) 
PUNCH  6»E(2)»E(3)»E(4),E(5)»E(6),E(7),E(8) 

"PUNCH  7 » F  mTFT2T TF (  3  ) TF T4)TFT5TTn~6TTFT7T  

PUNCH  9»P(1)»P(2)»P(3),P(4),P(5),P(6),P(7) 
"PUNCH  3 
IF    (JAC-7) 114, 111,93 
111  PUNCH  3 

PUNCH  3   

PUNCH  10 
PUNCH  11^ 
PUNCH  3 
PUNCH  12 
PUNCH  13 
PUNCH  3 

~~NDEC=T  

SUME=E( 1 ) 
SUMA=0. 
SUMP=0. 

DO   112  J=1,JAC   

DDK  =  D ( J ) /10. 
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Figure  5* — Continued 


112 


101 


93 


104 


105 


106 


107 


108 


109 


110 


114 


113 


DS=AREA*PK/100. 

DO   112  NDEC=1»10 

5UME=5UME+DDK 

SUMA=SUMA+DS 

SUMP=SUMP+PK" 

PUNCH  46  ♦SUME ,SUMAtSUMP 

IF    ( KACC-NUM) 100» 113 ♦ 113 

READ  45 »E<  1  )  *E( 2  )  »E( 3)  »E( 4) »E( 

READ  45»R( 1) ,R<2) »R(3) ,R(4) ,R{ 

READ  45»E(8)  »E(9)  *E(  10)  »E( 11)  ♦ 

READ  45  ♦  R  (  8  )  ,R(9)  »RTToT7RTTTJ  ♦ 

GO  TO  99 

PUNCH  3 

JA=JAC-7 


5 ) »E(6) »E<7) 
5 ) ,R(6  }  »R(7  ) 
E(12)  ,E(13)  tE(  14) 
RTY2T7R  TT37TR1  14) 


GO  TO 
PUNCH 
PUNCH 
PUNCH 
PUNCH 
GO  TO 
PUNCH 
PUNCH 
PUNCH 
PUNCH 
GO  TO 

PUNCH 
PUNCH 
PUNCH 
PUNCH 
GO  TO 
PUNCH 
PUNCH 
PUNCH 
PUNCH 
GO  TO 
PUNCH 
PUNCH 
PUNCH 
PUNCH 
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